Protection of center-spin coherence by a dynamically polarized nuclear spin core 
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Understanding fully the dynamics of coupled electron-nuclear spin systems, which are important 
for the development of long-lived qubits based on solid-state systems, remains a challenge. We 
show that in a singly charged semiconductor quantum dot with inhomogeneous hyperfine coupling, 
the nuclear spins relatively strongly coupled to the electron spin form a polarized core during the 
dynamical polarization process. The polarized core provides a protection effect against the electron 
spin relaxation, reducing the decay rate by a factor of Ni, the number of the nuclear spins in the 
polarized core, at a relatively small total polarization. This protection effect may occur in quantum 
dots and solid-state spin systems defect centers, such as NV centers in diamonds, and could be 
harnessed to fabricate in a relatively simple way long-lived qubits and quantum memories. 
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Introduction. — Controlling and extending the coher- 
ence time of a qubit lie at the heart of spintronics, quan- 
tum computation, and quantum information processing. 
Decoherence occurs inevitably because of the interaction 
of the system with its environment, which eventually 
makes the system behave classically. To counter with 
the decoherence, many proposals have been developed, 
including dynamical decoupling 1 -, decoherence-free sub- 
space^ 3 -, and environmental state preparation^—. 

Isolated electron spins in solids are promising infor- 
mation processors in spintronics and quantum comput- 
ing^, due to their long relaxation time in high mag- 
netic fields as demonstrated in systems of quantum dots 
(QDs)£. In low magnetic fields, however, rapid relax- 
ation may be caused by flips with surrounding nuclear 
spins, typically in a timescale of 10 ns in a few mT mag- 
netic field and at ^100 mK low temperature 1 ^—, which 
is not significantly longer than the nanosecond operation 
clock in electrical contro l 14 ! 15 . Thus extension of elec- 
tron spin coherence time is highly desired^. A partic- 
ularly promising scheme is to prepare the nuclear spins 
in low-fluctuation states^iiZ^— , which also supplements 
to the dynamical decoupling method by alleviating the 
requirements of pulse control. 

Nuclear spin state preparation has been investigated 
both theoretically and experimentally. Previous inves- 
tigations show that a nearly 100% nuclear spin po- 
larization, which is difficult to realize in experiments, 
is required in order to extend significantly the coher- 
ence time^ i 13 ' 23 : 24 , if a strong magnetic field is employed 
to uniformly polarize the nuclear spins. By contrast, 
the electron spin coherence time can be extended 10- 
100 times by dynamically polarizing nuclear spins us- 
ing repeated injection of polarized electron spins at a 
rather low total nuclear polarization of ~ l % 19 i 20 i 25 " 2E, 
But fully understanding the microscopic mechanisms of 
the dynamic nuclear polarization (DNP) of the coupled 
electron-nuclear spin systems remains a major challenge. 



Aiming at understanding the evolution of the coupled 
electron-nuclear spin system, we present in this paper a 
microscopic picture of the formation of a polarized core 
of nuclear spins during the DNP process by using a core- 
skirt two-region model [see Fig. QJa)]. We focus on the 
protection of the electron spin coherence by the polarized 
nuclear core in QDs, so as to provide a detailed and clear 
physical picture for the DNP effect. Our results show 
that the electron spin relaxation time can be extended 
substantially with a relatively small nuclear spin polar- 
ization. In a QD with inhomogeneous electron density, 
the strongly coupled nuclear spins (core spins) are eas- 
ier to be polarized during DNP than the weakly coupled 
ones (skirt spins). Once the N% core spins are polarized, 
they form a compound with the electron spin, which is 
decohered as a whole by the skirt spins. When the com- 
pound exchanges one quantum of moment with the skirt 
spins, the center spin moment is changed approximately 
l/(iVi + 1) fraction of a quantum. So the center spin 
relaxation is roughly Ni + 1 times slower (see Fig. 

Formation of polarized core during DNP. — We con- 
sider a coupled electron and nuclear spin system with 
a Hamiltonian (Gaudin spin-star model) 



V. = uj S z 



N 

£ 
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A k I k , 



(1) 



where ujq is the Zeeman energy of the central electron 
spin S {S = 1/2) and I k (I k = 1/2, k = 1, 2, • • • , N) is 
the fcth nuclear spi n 11 ' 12 : 30 . We set H = 1 for simplicity. 
Specifically, for an electron in a QD, A k is the Fermi con- 
tact hyperfine coupling constant, which is proportional to 
^(xfc)! 2 , the electron density at the fcth nucleus. Due to 
the inhomogeneity of \<ft\ 2 , A k is non- uniform in general. 
Such an inhomogeneity of A k is important to the forma- 
tion of polarized nuclear spin core during DNP. 

To illustrate the mechanism of formation of the po- 
larized nuclear spin core, we first consider a two-region 
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FIG. 1: (Color online) (a) Diagram of inhomogeneous hyper- 
fine coupling in a quantum dot and the polarized nuclear spin 
core (spins inside the dashed circle), (b) Time dependence 
of the polarization ratio r in a single DNP cycle for different 
A2 — 0.01, 0.05, 0.25, 1 from top to bottom. Other parameters 
are u = 6, Ni = 8, Na = 16, A\ = 1. The time unit is l/bi 
with bi — \J N\A\ . (c) Polarization ratio r versus the num- 
ber of DNP cycles. Other parameters are the same as in (b). 
Vertical dashed lines denote the cycle number 30, 40, 50. (d) 
Polarization ratio r versus hyperfine interaction strength ratio 
A2/A1 at DNP cycle 30 (circles), 40 (crosses), 50 (triangles). 



model in which the nuclear spins consist of a core and a 
skirt region. The hyperfine interaction strength Ai of N\ 
core nuclear spins is stronger than the strength A2 of N% 
skirt nuclear spins [see Fig. Ufa)]. So the Hamiltonian in 
Eq. UJ can be written as 



n = wo5, + i4iS-Ii + i4 a S-I a . 



Here Ii(2) = Ss^a' ^Jfe ^ s the total spin of the nuclei in 
the core (skirt) region. This Hamiltonian conserves I\, 
I2, and the z-component of the total spin S z + l\ z + l2 Z , 
which enables an exact formulation of the DNP process. 

Below we show that in a finite external magnetic field 
the polarization of the core spins acquired during the 
DNP process is much larger than that of the skirt spins, 
provided that Ai ^> A2. We assume that the initial elec- 
tron spin state is spin-up and the initial nuclear spin state 
is maximally mixed, p = 2- (Nl+N ^\ f)(t I ® 1 with 1 
being a unit matrix of dimension 2 Nl+N2 . Given that 
the electron-spin-mediated indirect coupling between Ii 
and I2 is negligible in a magnetic field ojo ^> y/N\$A\ t 2, 
the nuclear spin polarization saturates at a value p\.2 — 
(hz,2z) /Ni^ oc (Ai^/i^o) 2 at l° n g times, according to 



(2) 



the perturbation theory. Thus the polarization ratio 
r = Vilvi is proportional to the square of the local hyper- 
fine coupling strength, i.e., r oc (A1/A2) 2 . As shown in 
Fig. [JJ exact numerical calculations according to Eq. @ 
agree with the perturbation theory results. To initiate 
a new DNP cycle, we reset the electron spin to the up 
state, i.e., set the system to | f)(t l®Trs[p(T)] where the 
partial trace is over the electron and T is the cycle period. 
After many DNP cycles, numerical results in Fig. [JJ show 
that the relation r oc (A1/A2) 2 holds even better due to 
the fact that the total nuclear spin polarization (thus the 
effective ojq) increases with the number of DNP cycles. 
Of course, if the number of DNP cycles goes to infinity 
(much greater than N) , the maximum polarization of the 
nuclear spins is in the order of 1/\N for uniform hyper- 
fine coupling (Ak is constant) and reaches the order of 
1 for inhomogeneous AkS. Notwithstanding this limiting 
situation, the relation pk oc A^ usually holds in realistic 
experiments since the number of DNP cycles is < i V 19 i 29 . 

Coherence protection effect of the polarized core. — 
From the physical picture of the two- region model, we de- 
duce that when the electron-mediated indirect coupling 
between nuclear spins is negligible, the relation pk oc A\ 
should hold for a spin-star model with general inhomo- 
geneous coupling. To verify this, we perform numerical 
simulations with a Gaussian distribution of the hyperfine 
coupling coefficients Ak ■ Such a distribution can be real- 
ized, e.g., in a QD with a harmonic trap potential. The 
results (not shown) confirm the conclusion. 

To investigate the effect of DNP on the electron spin re- 
laxation, we simulate the dynamics of a coupled electron- 
nuclear spin system with Gaussian distributed Ak [see 
Fig. Ufa)] • To take into account the effect of nuclear spin 
polarization after many DNP cycles, we assume that in 
the initial state each nuclear spin has a polarization of 
Pk(fl) = ta,nh((3Af) with an adjustable parameter (3 re- 
lated to the inverse spin temperature^. We have ne- 
glected the possible phase correlations between nuclear 
spins built during the DNP process, which should decay 
much faster than the non-equilibrium spin polarization. 
For (3A\ < 1, pi- a A^, which reproduces the polar- 
ization ratio due to DNP with cycle number < N . As (3 
approaches infinity, pk saturates at 1, which accords with 
the limiting cases of an infinite number of DNP cycles. 

We consider both small- and large-bath cases, with 
the number of nuclei N = 20 and 256, respectively. 
For TV = 20, an exact evaluation is obtained using the 
Chebyshev polynomial expansion of the evolution oper- 
ator U = exp(-itH)2^; For N = 256, the method is 
based on P-representation of the density matri x 32 ' 33 . For 
N = 20, we also compare the P-representation method 
with the exact solution, and the results [Fig. Hfc)] show 
that the P-representation method provides a good ap- 
proximation. 
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FIG. 2: (Color online) Normalized decay of the electron spin polarization at various nuclear spin polarizations which increase 
from bottom to top for (a) N = 20 and (b) N = 256. The dashed horizontal lines mark the half decay point. The (last if 
multiple) cross point of the dashed horizontal line with each decay curve gives the half decay time Ti/a- (c) Dependence of 
Ti/2 (left j/-axis) and the nuclear spin core polarization (right y-axis) on the nuclear spin polarization in the case of TV = 20. 
Coincidence of the formation of the polarized core and the jump of the relaxation time T\i% manifests the key role played by 
the core. Black solid line with squares — the results of Chebyshev polynomial expansion method, black crosses — the results 
of P-representation method, dashed lines — core polarization, dash-dotted lines — skirt polarization, (d) Same as (c) for 
TV = 256. Black solid line with circles — the results of P-representation method. The Gaussian width is a for N — 20 and 6.4a 
for TV = 256, respectively, with a the lattice constant. The Gaussian center is shifted to (0.1, 0.27). 



Figure [2] presents the relaxation of an electron spin 
from the initial state | f). We set wo = to rule out 
the external magnetic field effect. The typical timescale 
is defined as b^ 1 = We show the results 

up to about 10 2 b^ 1 since other mechanisms (such as 
dipolar nuclear spin interactions) need to be included at 
longer times. The decay of the electron spin polarization 
is normalized as 

_ (S z (t)) - (5,M) 

- (5,(0)) - (S z (oo)) [ > 

where (S z (t)) = Ti{S z p(t)} with p(t) the density matrix 
of the whole system. 

In Fig. [21 significant extension of the relaxation time 
is observed with increasing the nuclear spin polarization. 
For N — 20, the decay time is extended by about 40 times 
when the nuclear spin polarization is changed from to 
70%. While for N = 256, a similar extension requires less 
than 20% nuclear spin polarization. Such a large elonga- 
tion of the decay time with a relatively small polarization 
contrasts sharply to the case of thermally polarized nu- 
clear spins, where polarization as high as 90% shows no 



significant extension of the decay tim o 12 i . The contrast 
between the DNP and thermal case indicates that the 
inhomogeneity in the nuclear spin polarization is of key 
importance to the extension of the decay time. In Fig. [21 
we also plot the polarization of 4 central nuclear spins 
in the core region. We see that nearly full polarization 
of the core spins coincides with the abrupt increase of 
the electron spin relaxation time. Actually, the factor of 
extension is roughly the number of nuclei in the nearly 
fully polarized core region. This phenomena suggest that 
the formation of a polarized core plays a critical role in 
protecting the center spin coherence. 

To understand the "protection" effect, we resort 
to the previous two-region model. If the core spins 
are fully polarized, they form together with the cen- 
ter spin a compound of N± + 1 polarized spins at 
the state |^tl,^ti) := \S, = §) ® \l u = f ) . 
After a flip with a skirt spin, the state of the 
core compound will be changed approximately to 

IJVi + l iVi-l \ / 2JVi+l I c _ 1\ (O, I T - ih-l\ , 

\^T~ '~2~ / — V 2(iVi + l) P z _ 2/ ^ \ llz ~ ~~2~/ + 

—r^ \S Z = -s)® \h z = ^) ■ Due to the hybridiza- 

. h(N, 4-11 I 2/ I lz 2 / J 
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FIG. 3: (Color online) Nuclear spin core protection effect, (a) 
Solid line denotes the electron spin polarization and dashed 
line denotes core spin polarization for the N — 20 case at 
68% total nuclear polarization [the largest polarization data 
in Fig.[2jc)]. (b) Two region model results. Solid lines denote 
the electron spin polarization for core spin number being Ni = 
(black line with circles), 2 (blue line with crosses), and 4 (red 
line). Dashed lines denote the core nuclear spin polarization 
for Ni — 2 (blue dashed line with triangles) and Ni = 4 (red 
dashed line). Other parameters are A\ = 1,A 2 = 0.1, N 2 ~ 
40, co = 0. The time unit is 1 /b 2 with b 2 = ^JN 2 A\. 

tion of the polarized core and the center spin, the electron 
spin is flipped with a probability of about l/(Nx + 1). 
Thus the relaxation time is extended by a factor of 
[N x + 1). 

To verify this picture, we examine the time-dependence 
of the nuclear core polarizations during the electron spin 
relaxation process. Figure[3Ja), which is numerically cal- 
culated for N = 20 with Gaussian distributed coupling, 
shows that the core polarization decays in accompany 
with the electron spin relaxation. The results for the 
2-region model shown in Fig. [3jt>) are similar. For dif- 
ferent number of core spins, the electron spin decay time 
increases linearly with the number of core spins N\ . All 



these features are consistent with the compound-spin pic- 
ture. 

Discussion and conclusion. — In a real QD system, the 
number of nuclear spins could reach as many as several 
millions^. According to the two-region model, 100 times 
extension of the coherence time requires a polarized core 
of about 100 nuclei, which corresponds to a total po- 
larization of ~ 0.01%. Considering the inhomogeneity 
of the hyperfine-interaction strength, we find it is very 
likely to extend ~ 100 times the electron-spin-coherence 
time with ~ 1% total polarization^. Full calculation of 
a double QD system will be presented in a future work. 

In conclusion, we show that in an electron-nuclei spin 
system with inhomogeneous coupling, the DNP can lead 
to the formation of a highly polarized nuclear spin core, 
which in return suppresses the electron spin relaxation 
with a relatively low degree of total nuclear spin polar- 
ization. Such effect may be observed in quantum dots, 
solid-state defect centers, and solid-state biomolecular 
NMR experiments^. In addition to the protection of 
the electron spin coherence, the polarized nuclear spin 
core is ready to be utilized to realize long-lived quantum 
memory based on imprinting and readout of electron spin 
state onto nuclear spins, which has a relaxation time on 
the scale of ten seconds— as proposed by Taylor and 
coworkers—. We anticipate that our results will also be 
of relevance in the full understanding of electron medi- 
ated nuclear spin diffusion during the DNP process of 
double QD systems^. 
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